Mesoscopic Physics of Granular Flows 
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We present a description of granular dynamics based on 
the idea of differentiation between fluid and solid components. 
First, we construct a model of completely fluidized phase. 
Then we discuss a shear surface motion on the boundary of 
the bulk solid phase, induced by a moving wall. Our results in- 
clude the thickness dependence of the velocity of a developed 
avalanche, density and velocity profiles in granular Couette 
flow experiment. 

PACS number: 45.70.M, 83.70.F, 05.20.J 

Dynamics of granular materials have recently attracted 
substantial attention of physicists M. Despite the nu- 
merous experimental and numerical results, there is still 
very little established theoretical framework for descrip- 
tion of granular flows. One of the dominant approaches 
is based on the principles of the kinetic theory of liq- 
uids S. Within this paradigm, it is believed that the 
dissipation in a granular system does not alter the basic 
structure of conventional hydrodynamics: the motion is 
expressed in terms of local velocity, density and (granu- 
lar) temperature. The dissipation can be accounted for 
by adding an appropriate dissipative term to the energy 
and momentum balance equation. Below, we argue that 
such an approach can hardly be valid for many granular 
systems, and propose an alternative approach to granular 
dynamics. 

Implicitly, kinetic approach assumes that energy and 
momentum in the system is transfered through the colli- 
sions between particles. In the case of hard-core interac- 
tions, these collision events take essentially no time com- 
pared to the free flight of the grains. However, this pic- 
ture does not survive in a general case of a dense dissipa- 
tive granular system. It is confronted by the numerically- 
discovered phenomenon known as inelastic collapse Q: if 
we start with an arbitrary distribution of grain velocities, 
the system of inelastic grains undergoes infinite number 
of collisions within a finite time and a large part of it 
virtually stops and creates dense clusters. 

If the system were sheared, the collapsed structures 
would provide a skeleton for the solid-like momentum 
transfer trough inter-grain mechanical contacts. This 
mechanism is supported by direct observation of force- 
carrying mesoscopic structures in the sheared granular 
system {H|. These string-like structures closely resem- 
ble the force chains responsible for stress transmission 
through static granular material and they are obvi- 
ously inconsistent with the collisional picture for the mo- 



mentum transfer. Our argument does not rule out the 
applicability of the kinetic theory to, for example, suffi- 
ciently dilute granular systems with weak dissipation, or 
to granular materials subjected to the bulk pumping of 
energy. 

In this letter, we are interested in physics of steady 
granular flow in the limit of high dissipation. As it was 
pointed out, in this case the fluid-like motion appears to 
be consistent with solid-like mechanism for stress trans- 
mission. The fundamental reason for this coexistence is 
that even if the contact network becomes percolating and 
can transfer the stress across the system, its connectivity 
may still be insufficient for mechanical stability (i.e. it 
may be below the rigidity percolation threshold). Fol- 
lowing this observation, we can distinguish between fluid 
and solid states of granular matter, based on its local 
connectivity. If the connectivity is high, the granular 
matter will be "jammed": a bulk of particles is only able 
to move as a solid, no motion of granules with respect to 
each other is possible. If, on the other hand, the connec- 
tivity of the network is low, grains can move with respect 
to each other. 

Thus, we can introduce a local "order parameter" 
as a local fraction of particles belonging to a given solid 
cluster. The motion of particle belonging to the same 
solid cluster is coherent, which makes the (\1/ = l)-phase 
similar to more conventional correlated states known in 
condensed matter physics, such as superconducting or 
superfluid condensates. In the following, it will be more 
convenient for us to use mobility parameter tf> = 1 — ^, 
which is the local fraction of the fluid component. Similar 
ideas of characterization of the granular system with the 
amount of fluidized component have been previously suc- 
cessfully used within BCRE model of surface avalanches 
0. The distinctive feature of our approach is that we 
are interested in the behavior of the system on the length 
scale of several beads, rather than in its macroscopic phe- 
nomenology. 

We limit ourselves to the generic case of planar ge- 
ometry, in which all principle variables, such as time- 
averaged velocity or the order parameter, vary only in 
the direction z normal to that of motion (x is directed 
along the flow). We shall consider two examples: flow 
of granules down an inclined plane and a shear flow of 
granular matter in Couette-type experiment. In the first 
case, the motion occurs under a given local stress, so 
the principle equations will include forces. In the second 
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case, the velocity on the inner wall is imposed, which is a 
kinematic condition, so the equations are also kinematic. 

Completely fluidized case: <f> = 1. We consider the 
motion of the granular material down an inclined plane, 
whose tilt is large enough for the system to be completely 
fluidized. In this state, the force chains transmitting an 
external stress are unstable with respect to buckling. In 
other words, the difference between a static force chain 
and one in the flow is that the latter has "weak bonds" 
to be broken due to external forcing (see Figure |l],a). 



"solid" 



its a zz and o~ zx components (given that o~ zy = 0). Now, 
we can rewrite the earlier result in the following form: 




b) 
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FIG. 1. Illustration to our discussion of flow mechanisms 
in "fluid" (a) and "solid" (b) phases. 

The network rearrangement associated with any of the 
buckling modes takes a finite time, which can be esti- 
mated by noticing that if / is the typical force transmit- 
ted by a force chain, the unbalanced torque at a particle 
of diameter d should be of order of fd. Over the buckling 
time, this torque should rotate the particle (of mass m) 
by an angle of order of unity, until a new contact is cre- 
ated. If Tb is the buckling time, this gives md 2 /t 2 ~ fd 



(1) 



The above estimate yields the time scale over which a 
single "weak" bond is destroyed. In a steady flow there 
is certain equilibrium fraction of the "weak" bonds (their 
number is equal to the number of independent "zero 
mode" deformations permitted by the contact network). 
This fraction can be expressed as a = d/£, where £ is 
the typical length of a stable chain segment between two 
unstable ("weak") points. Within our model, these sta- 
ble segments move like a solid, while the shear is con- 
centrated near the "weak" bonds. Therefore, £ can be 
interpreted as coherence length associated with the order 
parameter parameter ^F The shear rate d z v, is the time 
scale over which any particle will change its neighbors. 
We can estimate this rate as a jr. Now the shear rate can 
be related to the local stress by: d z v ~ a^ad/m, where 
a is the typical value of the stress. In the spirit of recent 
ideas on the constitutive law in granular solids ^,|), we 
assume that the stress tensor is completely determined by 



\d z v\ d z v 
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This equation of velocity has to be accompanied by an 
equation for the parameter a or, equivalently, coherence 
length £. It is a natural assumption that £ is a decreas- 
ing function of the ratio of the tangential and normal 
stresses. In fact, as this stress ratio decreases, the num- 
ber of free modes goes to zero, and £ should diverge at 
certain critical point, o zx j a zz — /i - After this jamming 
transition, the granular system becomes solid, i.e. a = 0, 
and £ = oo. 

If we start with the solid phase at the jamming point (mq 
and gradually increase the ratio a zx /a zz , e.g. by chang- 
ing the slope of the free interface with respect to the 
gravity direction, the solid will remain stable until an- 
other critical point, fj,\ > y«o. This means that £ is a 
two-branch function of the stress ratio: between /l*o and 
fii, £ = co for = 0, and £ is finite for 0=1. This 
non-uniqueness and the associated hysteresis, result in 
the stick-slip response to certain class of external driv- 
ing, observed experimentally [fiof . 

We now apply (Q) to the granular flow induced by an 
external bulk force (e. g. gravitational), tilted by an- 
gle 9 with respect to the z— direction. The breakage of 
bonds is dependent on the stress ratio <j zx /<j ZZi which re- 
mains constant everywhere and equals to tan 8. Let the 
height of the moving granular layer be h, with z = cor- 
responding to its bottom and z = h being its top. The 
component a zx of the stress is essentially given by the 
"hydrostatic" pressure: a zx — pg(h — z) sin#. Thus, (||) 
results in the following velocity field : 



v(z) 



Vgsin(6>) 
£(tan0) 



(h 3 ' 2 -(h- z) 3 ' 2 ) 



(3) 



We conclude that the typical velocity of a thick developed 
avalanche scales as h 3 ! 2 with its height, and that the 
dependencies of the average flow velocity on the layer 
height collapse onto a single master curve 



(4) 



for a variety of physical parameters (i.e., slope, friction, 
etc). The scaling law v ~ h 3 / 2 have indeed been ob- 
served experimentally by Azanza et al Moreover, 
Pouliquen |12| in recent study has demonstrated the va- 
lidity of above master relationship, Eq(Q). In the experi- 
ment, the fundamental length scale £ reveals itself as the 
minimal layer thickness at which the material flows for 
a given slope, which is completely consistent with our 
interpretation of £. The experiments also support our 
observation that £ diverges as the slope approaches the 
critical value, tan -1 
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Partially fluidized case (0 ~ 0). Our further dis- 
cussion is devoted to the shear dynamics of the solid 
phase perturbed by a moving wall, as in Couette flow 
experiment. Since granular system tends to create col- 
lapsed solid phase, the shear motion normally occurs only 
within a mesoscopic layer (several bead diameters wide) 
near the moving interface. Attempts to construct contin- 
uous description may be inadequate on this length scale. 
That is why we make z a discreet variable, e.g. imagine 
that the beads near the surface are organized in single- 
bead-wide layers (which in reality are reasonably well- 
defined). The layer at z = belongs to the moving wall, 
and has velocity vq; in the limit z = oo all beads belong 
to a single coherent solid cluster and have zero velocity. 

Consider n-th layer. The typical rate with which a 
particle belonging to it is being hit by a particle at the 
adjoint layer, n — 1, is v n —\jd. Any time when this hap- 
pens, there is a chance that the bead starts moving (gets 
"excited" ). The particle will move by a certain distance 5 
(its free path) until being stopped by another particle at 
the same layer, as shown on Figure 0(b) (one can easily 
modify the model to account for the possibility of "sec- 
ondary excitations" to be born). As a result, the average 
velocity of the n — th layer is 



a(v n -i + v n+ i) 



(5) 



Here a is the coupling coefficient which measures the 
probability for a particle to be moved after being hit by 
a neighbor. We made sure that v n -\ and v n +± appear at 
the above equation in a symmetric manner. The asym- 
metry comes in boundary conditions: once the position 
of the source of the excitations (the moving wall) is spec- 
ified, the velocities decay very rapidly with n and one of 
the velocities becomes negligible compared to the other. 

The obtained equation does not give a complete de- 
scription of penetration of the motion inside the solid 
phase, because nothing was said about the typical inter- 
particle gap (S) n . It would be reasonable to assume that 
these gaps make a leading contribution to the free vol- 
ume of the system, and therefore, (6) can be related to 
the deviation of the particle density p n from its random- 
closcd-packing value p cp : 



Pep 



-d. 



Pep 
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Eqs. (H|) and (0) provide and opportunity to relate den- 
sity and velocity profiles obtained experimentally. Note 
that constant (8) n would result in exponential depth de- 
pendence of the velocity. It is more natural to expect 
the available free volume to decrease with depth n, thus 
resulting in a stronger than exponential velocity decay. 

We present one possible way of modeling the depth 
dependence of S. We shall call this treatment the Quan- 
tized Free Volume approach. It is based on the observa- 
tion that after a particle jumps and gets stopped by its 



neighbors, the free volume available to it becomes zero. 
As a result, we expect the free path parameter 5 to be 
either zero (for particles in contact with their front neigh- 
bors) or of order of d (for movable particles), as shown 
on Figure |l|(b) . In the spirit of our earlier discussion, we 
interpret those movable particle as fluid component, and 
identify their local fraction with <f>. If 8„ ~ d is a mean 
free path for the fluid sub-system, then 



(S) n = 5*4>n 



(7) 



In order to find a missing equation for 5, or <fi, we 
consider an exchange of free volume between adjoint lay- 
ers. In the limit of <fi <C 1, our assumption that the 
free volume is quantized, has an important implication. 
Since inter-layer free volume exchange may occur only 
near movable particles (free volume quanta), its flux is 
linear with their concentration. The natural time scale 
of the volume exchange between n-th and n + 1-th lay- 
ers is given by v n —i/d (the rate of change of the near- 
est neighborhood). Hence, the general expression for the 
corresponding free volume flux is 



7 _ Vn ~ 1 i A 

Jn,n+1 — — J — {0~+9n 



(8) 



V n+ 2 



Here u n -|_2-term is added in order to make the structure 
of the model consistent with the inversion symmetry. In 
reality, after the boundary conditions are chosen, one of 
the two terms can be neglected, since velocity decays 
very rapidly with n. In the case of shear flow, we can 
neglect u n +2 Wn-i- I n t ne steady state, the flux given 
by (||) should vanish. The condition that J n ,n+i = 
provides the desired closure for (^-0). The result is an 
exponentially-decaying 0-profile: 



t> exp 



("5 



(9) 



Here A = l/(log(cr_/cr+)) is the penetration depth of the 
fluid component into the solid bulk. We assume a + < 
<7_, since otherwise the solution (^|) is clearly unphysical. 
Parameter (f>o ~ 1 depends on the detail of interactions 
between the moving wall an the first layer. Given the 
solution (^|), we can determine velocity field from (||) 



Vq exp 



1 

2A 



n+- 



X 

A* 



(10) 



Here A* = l/(\og(d/8^(j) a)) is another length scale in 
the problem. Its variability from one system to another 
is expected mostly through the coupling coefficient a. In 
the limit of very weak inter-layer coupling a (which corre- 
sponds to very smooth particles), A* may become much 
less than A, which would result in exponential velocity 
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profile: v n ~ exp(— n/A*). If the coupling is strong, both 
lengths A and A* are expected to be of the same order, 
thus making the velocity profile Gaussian shaped. The 
analytic and numerical results of our model are presented 
on Figure both for strong and weak coupling regimes. 
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FIG. 2. Velocity profiles (circles) and 0-field (diamonds) 
in partially-fluidized system. White symbols correspond to 
the weak coupling regime (A* = 0.4), and black ones corre- 
spond to the strong coupling (A* = 4.5). One can clearly 
see the crossover from exponential to Gaussian velocity pro- 
file. Curves and data points represent analytic and numeric 
results, respectively. A = 2 for all curves. 

Existing experiments do show super-exponential decay 
of the velocity profiles, in some systems it is remarkably 
close to Gaussian (5]Q| . Our model suggests the transition 
from exponential to Gaussian velocity profile depends on 
the strength of inter-layer coupling a. Note that this 
result differs from the conclusions of recent Debregeas- 
Josserand model ]l3| , which attributes the change in ve- 
locity profile to the effect of dimensionality. 

In conclusion, we have presented a description of gran- 
ular flows based on the notion that dissipation results 
in a strong correlation of particles' motion. An extreme 
case of such correlated motion is a coherent solid phase 
(4> = 0). We have constructed a model for the partial sur- 
face fluidization of this phase near a moving wall. The 
basic structure of our model employs some concepts sim- 
ilar to those of the conventional condensed matter the- 
ory: the motion is associated with discreet excitations 
(movable particles) "penetrating" into the coherent solid 
phase within a mesoscopic surface layer. We made a 
quantitative prediction relating velocity and density pro- 
files (p[]lO|), which can be checked experimentally. One 
more check of our assumptions is possible by studying the 
velocity distribution function with high temporal resolu- 
tion. If our approach is valid, one should be able to iden- 
tify the condensed and fluidized sub-systems with such 
kind of experiment. In the case of completely fuidized 
state, 4> — 1 , the collective behavior reveals itself through 



existence of short-living force chains. Once again, it is in- 
structive to focus on the mesososcopic scales, set by the 
typical length of a stable force chain segment (coherence 
length £). Our approach gives several important results, 
which include v(h) scaling law for thick avalanches, as 
well as velocity and density profiles in granular Couette- 
flow experiment. Our conclusions are in a good agree- 
ment with existing experiment, and some of them are 
yet to be tested. 

ACKNOWLEDGEMENTS 

We thank Georges Debregeas, Cristophe Josserand, 
Dan Mueth, Sidney Nagel, Heinrich Jaeger, Leo P. 
Kadanoff, Tomas A. Witten for fruitful discussions. This 
work was supported in part by MRSEC Program of Na- 
tional Science Foundation under Award Number NSF 
DMR-9808595. 



[1] H. Jaeger, S. Nagel and R. Behringer, Rev. Mod. Phys. 

68, 1259 (1996). 
[2] S. B. Savage and M. Sayed, J. Fluid Mech. 142, 391 

(1984). 

[3] I. Goldrisch and G. Zanetti, Phys. Rev. lett. 70, 1619 
(1993). 

[4] C. Veje, D. Howell, R. Behringer, S. Schollmann, S. Lud- 
ing and H. Herrmann, in Physics of Dry Granular Media, 
NATO ASI Series E350, 237 (1998); C. T. Veje, D. W. 
Howell, and R. P. Behringer, Phys. Rev. E 59, 739 (1999). 

[5] D. M. Mueth, G. F. Debregeas, G. Karczmar, H. M. 
Jaeger and S.R. Nagel, to be published. 

[6] G. W. Baxter, Powders and Grains 97, (R. P. Behringer 
and J. T. Jenkins, eds), Rotterdam, Balkema (1997); C. 
Liu, S. Nagel, D. Schecter, S. Coppersmith, S. Majum- 
dar, O. Narayan and T. Witten, Science 269, 513 (1995). 

[7] J.-P. Bouchaud, M. Cates, R. Prakash, and S. F. Ed- 
wards, J. Phys. (France) 4, 1383 (1994); J.-P. Bouchaud, 
M. E. Cates, J. R. Prakash, and S. F. Edwards, Phys. 
Rev. Lett.74, 1982 (1995). 

[8] J. P. Wittmer, P. Claudin, M. E. Cates and J. P. 
Bouchaud, Nature 382, 336 (1996); J. P. Wittmer, M. 
E. Cates and P. Claudin, J. Phys. I France 7, 39 (1997). 

[9] A. V. Tkachenko and T. A. Witten, Phys. Rev. E.60, 687 
(1999). 

[10] S. Nasuno, A. Kudrolli, A. Pak, and J. P. Gollub, Phys. 

Rev. E. 58, 2161 (1998). 
[11] E. Azanza, F. Chevoir, and P. Moucheront, Powders and 

Grains 97, (R. P. Behringer and J. T. Jenkins, eds.), 

Rotterdam, Balkema (1997). 
[12] O. Pouliquen, preprint. 

[13] G. Debregeas and C. Josserand, preprint |cond- 
mat/9901336[ 



4 



